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Abstract. The ISOPHOT instrument aboard ISO has been 
used to observe extended FIR emission of six Abell clus¬ 
ters. Strip scanning measurements with crossing position an¬ 
gles centered on the clusters were cari'ied out at 120/rm and 
180 pm. The raw profiles of the Ii 20 Atm/Ii 80 /im surface bright¬ 
ness ratio including zodiacal light show a bump towards Abell 
1656 (Coma), dips towards Abell 262 and Abell 2670, and are 
without clear structure towards Abell 400, Abell 496, and Abell 
4038. After subtraction of the zodiacal light, the bump towards 
Abell 1656 is still present, while the dips towards Abell 262 and 
Abell 2670 are no longer noticable. This indicates a localized 
excess of emitting material outside the Galaxy towards Abell 
1656 with properties different from the galactic foreground cir¬ 
rus, while the behavior in Abell 262 and Abell 2670 can be 
reconciled with galactic citTus structures localized on the line- 
of-sight to these clusters. The excess of « 0.2 MJy/sr seen 
at 120 pm towards Abell 1656 (Coma) is interpreted as being 
due to thermal emission from intracluster dust distributed in 
the hot X-ray emitting intracluster medium. The integrated ex¬ 
cess flux within the central region of 10' - 15' diameter is w 
2.8 Jy. Since the dust temperature is poorly constrained, only a 
rough estimate of the associated dust mass of Md ~ 10^ Mq 
can be derived. The associated visual extinction is negligible 
(Ay ^ O.lmag) and much smaller than claimed from opti¬ 
cal observations. No evidence is found for intracluster dust in 
the other five clusters observed. The absence of any signature 
for intracluster dust in five clusters and the rather low inferi'ed 
dust mass in Abell 1656 indicates that intracluster dust is likely 
not responsible for the excess X-ray absorption seen in cool¬ 
ing flow clusters. These observations thereby represent a fur¬ 
ther unsuccessful attempt in detecting the presumed final stage 
of the cooling flow material, in accord with quite a number of 
previous studies in other wavelengths regions. Finally, the ob¬ 
served dimming of the high-redshift supernovae is unlikely be 
attributable to extinction caused by dust in the intracluster or 
even a presumed intercluster medium. 

Key words: galaxies: clusters: general - galaxies: clusters: in¬ 
dividual: Abell 262, Abell 1656, Abell 2670 - intergalactic 
medium - infrared: general 
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1. Introduction 


The analysis of the X-ray and UV emission of the hot (Te « 1 
- 10 keV) el ectron gas trap ped in the potential well of galaxy 
clusters (e.g. Sarazin 1986| ) has revealed an increasingly com¬ 
plex picture of the intracluster medium (ICM). High resolu¬ 
tion X-ray images reveal the pr esence of small-scale structures 
of different temperatures (e.g. Donnelly et al. 1999, and ref¬ 


erences therein). High central densities, lower temperatures, 
and cooling times significantly shorter than a Hubble time in 
the centers of many clusters have been interpreted as evidence 
for cooling flows (e.g. Fabian 1994; White et al. 1997; Allen 
2000 ), where the cooling gas drops out of thermal equilibrium, 
possibly leading to the formation of dense clouds and even¬ 
tually to stars. X-ray spectra and image decomposition indi¬ 
cate an extended excess absorption of soft X-rays with an HI 
column density of « 10^^cm“^ above the galactic foreground 


(e.g. White et al. 1991; 

Allen & Fabian 1997; 

Sarazin 1997 

Allen 200CD. However, 

krabadjis & Bregman (2000) And sig- 


niflcantly lower or even no additional absorbing column den¬ 
sities in a number of other clusters. A spatially extended EUV 


emission of non-thermal origin (Enfllin et al. 1999, and ref¬ 


erences therein) above the expected flux from the hot X-ray 
gas in at least the Co ma and possibly o th er clusters seems to 
be fairly established (Lieu et al. 1996a ,^; Mittaz et al. 1998 ; 


Lieu et al. 1999 ; Bowyer et al. 1999 ), although an accompa¬ 
nying lo w energy X-ray excess emissi on is uncertain (Dixon et 
al. 1996; Arabadjis & Bregman 1999 ). Overall, these observa¬ 
tions indicate a striking departure from a nearly homogeneous 
intracluster gas towards a spatially inhomogeneous multi-phase 
ICM (Bonamente et al. 2001; ^ieu et al. 200C ). 

Despite quite some observational effort, neither the final 
stage of the cooling intracluster gas nor the material responsi¬ 
ble for the e xcess absorption has been detected at other wave¬ 
lengths (e.g. Sarazin 1997 ; ji^aor 1997| ; |Koekemoeret al. 1998| ; 
Miller et al. 1999; Allen 2000, and references therein). It has 
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been suggested that the absorbing material is present as either 
very cold molecular clouds ( |^erland et al. 1994 ) or dense cold 
clouds where the shielding allows the formation of dust grains 
onto which_alWolutilesJiaveftozen^ames_e£aLJ^9^Fabian 
et al. 1994). 

The measured line emission of heavy elements (e.g. 
Mushotzky et al. 1996 ; Fukazawa et al. 1998 ) clearly indicates 
that the ICM contains a processed component, likely the result 
of galactic winds, stripping of the interstellar medium (ISM) 
of cluster members by ram pressure, or tidally removed ISM 
during merging with other galaxy groups. These processes also 
add dust to the ICM, which in turn is a possible candidate for 
the excess X-ray absorber. The first direct evidence for dust in 
the ICM is an observed oxygen K edge in an X-ray spectrum 
of the Perseus cluster attributed to neutral material, most likely 
dust grains condensed out of the ICM (Arnaud & Mushotzky 
1998). 

In equilibrium conditions, the predicted intracluster dust 
(ICD) temperatures are rather uncertain with temperatures in 
the range 10 - 20 K (Dwek et al. 1990; Loewenstein & Fabian 
1990; Braine_e£^lM^95^j_but£ossibl^jnuchJiigher_(0’Dea et 
al. 1994; |Voit & Donahue 1995| ). The thermal emission of the 
ICD is thus to be expected in the FIR. If not shielded in dense 
clouds, the high kinetic temperatures destroy dust grains with 
a typical size of ~ 0.1 /im by sputtering o ver time scales of a 
few 10® - 10® years ( Dwek et al. 199C ; Helens et al. 1994 ; 
Dwek et al. 1996 ). Very big (> 10/im) dust grains introduced 
to account for the observed FIR spectra of galaxies (Rowan- 
Robinson 1992) would have a significantly longer lifetime of 
the order a Hubble time. Possibly, the dust destruction is more 
efficient for smaller (< O.l/rm) grains, leaving preferentially 
behind the bigger grains, which moreover might have a flatter 
extinction curve than those observed in the Milky Way (Aguirre 
1999b). Furthermore, the temperatures in the central cooling 
flow regions of clusters might generally be lower, which, in 
addition to the shielding by clouds, would contribute to an in¬ 
creased grain lifetime. 

If dust is the absorbing material inferred from the X-ray 
observations, a significant fraction of the absorbed X-ray radia¬ 
tion of ^ 10^® erg is re-emitted in the FIR , and dust would 
then constitut e a major coolant of the ICM ( Bregman 1992 ; 
ISarazin 1997). The FIR emission from ICD would then also 
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be relevant for the correct interpretation of measurements of 
the Sunyaev-Zeldovich effect in the sub-mm wavelength range 


Mpc, Hg = 70 km s ^ Mpc ^, qg = 0.5) diameter of the clus¬ 


(Khersonskii & Voshchinnikov 1985; Lamarre et al. 1998). 


ter was observed (Stickel et al. 1998). In a follow-up study, 


Intracluster dust may reveal itself by optical depth effects 
in cosmolog ical investigations of distant objects ( Shaver 1987 ; 
Masci 1998 ), a significant fraction of which are viewed through 
clusters. Historically, this has be en the first attempt to find ev¬ 
idence for the existence of ICD (Zwicky 1957 1962). Despite 
quite some effort to search for the visual extinction of back¬ 
ground quasars and galaxies seen through foreground galaxy 
clusters, no common consent has been achieved whether prop¬ 
erties of background sources are actually affected by the pres¬ 
ence of dust in foreground clusters (Girardi et al. 1992; Fer¬ 
guson 1993; Maoz 1995, and references therein). Aguirre 
(1999a,^ and Simonsen & Hannestad (1999) have drawn at¬ 


tention to the possibility that the observed dimming of distant 
type la supernovae might be caused by intervening dust rather 
than resulting from a non-zero cosmological constant. 

From a comparison of Ly a and UV continuum fluxes in 10 
cooling flow clusters, |du (1992 ) found an average excess red¬ 
dening of Eb-v ~ 0.19 above the galactic foreground absorp¬ 
tion. This was attributed to ICD spread out over a large fraction 
of the cluster volume. Model calculations of the extended FIR 
emission from ICD were presented by pwek et al. (199C ) and 
Hu (1992), who concluded that the expected FIR flux was just 


ISOPHOT was used to observe five more galaxy clusters in a 
similar way, and the results, together with a re-analysis of the 
Coma data, are described in the following. 


2. Cluster Selection 

In the hrst intracluster dust observation with ISOPHOT (Stickel 
et al. 1998), the Coma cluster was selected because it is the 
cluster where early optical extinction studies had found ev¬ 
idence for ICD ( ^wicky 1962| ; Karachentsev & Lipovetskii 
1969), where diffuse optical emission from intracluster mate¬ 
rial was already known (Welch & Sastry 1971; Mattila 1977), 
and where a detailed theoretical study ( Dwek et al. 199^ 
showed that the predicted diffuse FIR emission was within 
reach for ISOPHOT. 

Primary candidates for ad ditional observation s were Abell 
262 and Abell 2670, since Wise et al. (1993 ) had found 
marginal evidence for extended FIR emission on IRAS ISSA 
plates. Because the Coma clu ster is known to underg o merging 
with smalle r galaxy groups (Colless & Dunn 1996| ; Vikhlinin 
et al. 1997; Burns et al. 1994b) the dust inferred from the dif¬ 


within reach of the available IRAS data and the then upcoming 
ISO mission. 

The FIR emission as a direct indicator for ICD has been 
searched for at 60 //m an d 100 /im on IRA S ISSA plates 
(Wheelock et al. 1994) by Wise et al. (1993). Removing the 
large scale galactic foreground cirrus emission, evidence for 
diffuse excess FIR emission of low statistical signihcance was 
found in the direction of several out of 56 clusters. Sub¬ 
millimeter observations of 11 cooling flow clusters by Annis 
& Jewitt (1993) did not detect any emission from dust near the 
cluster centers. With the more sensitive sub-millimeter bolome¬ 
ter array at the JCMT, ^dge et al. (1999 ) detected the dust 
emission of the central galaxies from 2 out of 5 cooling flow 
clusters. However, the latter two studies were insensitive to ex¬ 
tended emission of ICD on angular scales larger than a few 
arc minutes. Hansen et al. (200^ used ISOPHOT FIR maps to 
search for dust coincindent with the central galaxies of cooling 
flow clusters, but obtained only inconclusive results. 

Compared to FIR maps, strip-scanning measurements 
crossing the cluster provide a means to search a larger frac¬ 
tion of the cluster volume, and to reach the FIR foreground 
level far off the cluster centre. Furthermore, the separation of 
extended intracluster dust emission from the foreground zodi¬ 
acal and galactic cirrus emission can be accomplished by us¬ 
ing observations at two different wavelengths longward of the 
IRAS bands. This combined technique was employed in a pi¬ 
lot study to search for dust in the intracluster medium of the 
Coma cluster (Abell 1656, z = 0.0232) with ISOPHOT (Lemke 


fuse FIR emission was interpreted as being transferred rather 
recently to the ICM by tidal stripping of cluster galaxy ISM 
(Stickel et al. 1998). Extended FIR emission from the ICM 


thus might be found preferentially in recent mergers or young 
clusters, which makes the testable prediction that old clusters 
should show weak or no extended FIR emission. It was there¬ 
fore attempted to And more targets along the sequence of X-ray 
morphologies outlined by lones & Forman (1992) and quan¬ 
tified as a possible dynamical and evolutionary sequence by 
Buote & Tsai (1995, 1996). Particularly, the supposedly young 
clusters with irregular and the supposedly old virialized clus¬ 
ters with circular symmetric X-ray emission were considered. 

Cluster candidates from the list of Buote & Tsai (199^ 
were subject to constraints in the sky visibility by ISO and in 
cluster redshift. The latter was necessary because at too large 
redshifts only very few non-overlapping sky positions would 
have been possible to observe with the 3' wide ISOPHOT de¬ 
tector inside a characteristic cluster Abell radius, while the 
nearest clusters would have required prohibitively long integra¬ 
tion times to cover a significant fraction of the cluster diameter. 
Since a re-observation of the Coma cluster became impossible, 
Abell 4038 was chosen as a closely resembling cluster, par¬ 
ticularly wit h respect to the mul tipole moments of its X-ray 
morphology ( Buote & Tsai 1996 ). 

The observed Abell clusters are listed in Tab. which 
gives for each cluster (Col. 1) the coordinates, on which the 
ISOPHOT observations were centered (Col. 2, 3), its galactic 
coordinates (Col. 4, 5), the cluster redshift (Col. 6), the cluster 
richness class (Col. 7), and the cluster Abell radius according 


et al. 1996; Lemke & Klaas 1999), the photometer instrument to Rudnick & Owen (1977 ) (Col. 8). 


aboard ISO ( Kessler et al. 199^ 


Kessler 1999| ). An enhanced 
120/xm emission within the central region of « 10' (« 0.3 
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Table 1. Observed Abell clusters 


Cluster 

Q2000 

^2000 

1 

b 

Redshift 

Richness- 

I^Abell 




n 

n 


Class 

['] 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

Abell 262 

0l'*52"*46.4'* 

-f36°09'06" 

136.58 

-25.09 

0.0161 

0 

110 

Abell 400 

02'*57’"38.6® 

-b06°02'00" 

170.24 

-44.94 

0.0238 

1 

75 

Abell 496 

04'*33"*37.1'* 

-13°15'42" 

209.57 

-36.48 

0.0327 

1 

56 

Abell 1656 (Coma) 

12'*59"*35.7'* 

-f27°57'38" 

58.16 

-f88.01 

0.0232 

2 

75 

Abell 2670 

23'*54"*13.7'* 

-10°25'09" 

81.32 

-68.52 

0.0761 

3 

26 

Abell 4038 (Klemola44) 

23'*47’"41.9® 

-28°08'19" 

25.08 

-75.90 

0.0283 

2 

64 


Notes : 


■ coordinates (Col. 2, 3) of central sky position of ISO scans 

■ redshifts (Col. 6) and richness cl asses (Col. 7) taken from NED 

■ Abell radii (Col. 8) according to Rudnick & Owen (1977) 


Table 2. Observational Details 


Cluster 

Observing 

Date 

ISO 

Revolution 

# Steps 

Cluster 

Coverage“ 

Position Angle 

n 

On-target 

Time 

[s] 

Zodiacal Light at 
120 pm / 180 pm 
[MJy/sr] 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

Abell 262 

Jan 7, 1998 

784 

31 

0.82 

45, 135 

10250 

2.7/1.8 

Abell 400 

Feb 21, 1998 

829 

19 

0.76 

40, 130 

6530 

5.1 / 2.4 

Abell 496 

Mar 12, 1998 

848 

17 

0.91 

90 

2950 

1.9/0.9 

Abell 496 

Mar 19, 1998 

855 

17 

0.91 

180 

2950 


Abell 1656 (Coma) 

Jul 21, 1996 

247 

16 

0.64 

36, 82 

3480 

2.3/1.1 

Abell 2670 

Nov 21, 1997 

736 

17 

1.96 

19, 84 

5640 

3.8/1.6 

Abell 4038 (Klemola 44) 

Dec 24, 1997 

769 

19 

0.89 

120 

3260 

3.0/ 1.4 

Abell 4038 (Klemola 44) 

Dec 25, 1997 

771 

19 

0.89 

30 

3260 



Notes : 

“ scan length in units of the Abell radius 


3. Observations 


The observations were carried out with the C200 camera of 


ISOPHOT (Lemke et al. 1996; Lemke & Klaas 1999), a 2x2 


pixel array of stressed Ge:Ga with a pixel size of 89"4. For each 
cluster, linear scans across the cluster center at two different po¬ 
sition angles were obtained, each of which was observed with 
both the C_120 and C_180 filters (reference wavelength 120 fim 
and 180 fim, equivalent widths 47 /im and 72 /rm, respectively) 
at the same sky positions. The positional offset between subse¬ 
quent sky measurements was 3', which resulted in almost no 
detector overlap. The integration time at each sky position was 
« 50s for Abell 1656 and ss 80s for the other clusters, dur¬ 
ing which 6-16 integration ramps with 127 non-destructive 
readouts were obtained. 

Details of the observations are listed in Tab. ^ which gives 
for each cluster (Col. 1) the date (Col. 2) and ISO revolution 
number of the observation (Col. 3), the number of sky positions 
(Col. 4) and the fraction of the Abell radius covered by the 
scan (Col. 5) at each position angle (Col. 6) together with the 
total on-target integration time (Col. 7). Finally, Col. 8 gives 
the zodiacal light contribution for 120 pm and 180 pm for the 


observing date in Col. 2 taken from Kelsall et al. (1998 ) for 
all clusters except Abell 2670, which in turn was taken from 
the yearly averaged zodiacal light map given by Leinert et al. 
(1998). 


Physical properties of the observed clusters are collected in 
Tab. H which gives for each cluster (Col. 1) the overall X-ray 
morphology (Col. 2), the galactic foreground column density 
(Col. 3), the X-ray temperatures (Col. 4), the estimated cool¬ 
ing flow rate (Col. 5), the derived gas mass (Col. 6), and the 
X-ray luminosity (Col. 7). Where possible, several entries for 
each quantity together with the literature reference have been 
listed to give an indication for the spread in the determination 
of these values. The six observed clusters cover the range from 
relaxed elliptical (Abell 262) to irregular (Abell 400) X-ray 
morphologies, from cool (Abell 262, Abell 400) to hot (Abell 
1656) X-ray temperatures, from no (Abell 400, Abell 1656) to 
a significant (Abell 496) cooling flow rate, and about an order 
of magnitude in galactic Nhi column density. 
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Table 3. Physical Properties of Observed Clusters 


Cluster 

X-Ray 

NhI 

Temperature 

Cooling Flow 

GasMass 

LxBoI 


Morphology 

[10^°cm-2] 

[keV] 

[M0 yr"^] 

[IO^^Mq] 

[lO'^'^ergs-i] 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Abell 262 

elliptical [4] 

5.3 [1] 

2.4 [1] 

27 [3] 

2.5 [3] 

0.9 [1] 



10.5 [4] 

2.5 [2] 

9 [2] 

6.5 [2] 

0.5 [2] 




1.4 [4] 

47 [7] 


0.4 [4] 

Abell 400 

irregular [6,8] 

8.8 [1] 

2.5 [1] 

0 [2] 

11.1 [2] 

0.7 [1] 




2.2 [2] 

0 [10] 

15.0 [6] 

0.6 [2] 

Abell 496 

single 

4.4 [1] 

3.9 [1] 

112 [7] 

2.7 [2] 

5.8 [1] 


symmetric [8] 

12.0 [5] 

4.8 [2] 

138 [2] 

74.1 [2] 

7.9 [2] 




3.3 [5] 

95 [3] 



Abell 1656 

elliptical [11] 

0.9 [1] 

8.3 [1] 

0 [3] 

27.8 [3] 

12.5 [1] 

(Coma) 



7.3 [2] 

0 [2] 

75.6 [2] 

9.1 [2] 

Abell 2670 

double [8,9] 

2.7 [1] 

3.9 [1] 

41 [2] 

40.7 [2] 

3.8 [1] 



3.6 [9] 

3.7 [2] 

0 [9] 


3.5 [2] 




4.2 [9] 

18 [12] 



Abell 4038 
(Klemola 44) 

elliptical [8] 

1.5 [1] 

3.3 [1] 

87 [3] 

4.3 [3] 

2.2 [1] 


References : 

[1] David etal. (1993); [2] White et al. (1997); [3] Peres et al. (1998); [4] David et al. (1996); [5] MacKenzie et al. (1996); [6] Beers et 
al. (1992); [7] Edge et al. (1992); [8] Buote & Tsai (1996); [9] Hobbs & Willmore (1997); [10] Allen et al. (1997); [11] Vikhlinin et al. 
(1997); [12] Wise etal. (1993) 


4. Data Reduction 


To push the detection of the expected weak signals at these 
long wavelengths to the limits, and to get confidence in de¬ 
tected features, several different methods for the signal deriva¬ 
tion we re used in parallel. Already in the case of the Coma 
cluster ( Stickel et al. 1998 ), the pairwise differences of con¬ 
secutive ramp read-outs rather than signals from the full ramp 
fitting have proven to provide the most robust signals, particu¬ 
larly if only a few ramps at each sky position are available. To 
get rid of pairwise readout differences affected by cosmic ray 
hits, the robust outlier-insensitive myriad estimator was com¬ 
puted and 15% of the most deviant signals as measured by the 
absolute deviation were cut off. This outlier removal is simi¬ 
lar to a median absolute deviation trimming, but instead of the 
initial median, the sample myriad is used to determine the out¬ 
liers. The sample myriad value in turn is a robust estimator of 
the mode (most common value) of a distribution but does not 
require binning of the actual data set, and is easily computed 
by minimizing a particular cost function with a tuning constant 


set to a small value (for details see Kalluri & Arce 1998). 


Separate sets of final signals for each sky position were 
derived from this tail-trimmed pairwise distribution by com¬ 


puting the mean, median, myriad (Kalluri & Arce 1998), and 


the annealing M estimator (Li 1996). Additionally, the cumula 


tive trimmed distribution was fitted with the analytical function 
y = Co(l + exp((x — C 3 )/ci))‘^^), the first derivative of which 
resembles very closely a gaussian. In this case, the final signal 


value was derived from the maximum of the first derivative of 
the fitting function. 


Signals from the much smaller distribution of ramp slopes 
at each sky position, obtained by fitting first order polynomials 
to the deglitched readouts of each integration ramp, were de¬ 
rived for comparison. The initial ramp-fitting step was done 
within the ISOPHOT Interactive Analysis PIA[] Version 9.0 
(Gabriel et al. 1997). Since the number of ramps for one sky 
position was rather small, a subsequent efficient rejection of 
disturbed ramp slopes was not possible. Therefore, only the ro¬ 
bust myriad, the annealing M estimator, and the standard PIA 
averaging were applied to derive the final signals. 


Eventually, the signals derived by all different methods 
were corrected for signal dependence on ramp integration times 
to be consistent with calibration observations (Laureijs et al. 
2000), dark-current subtracted, and finally flux calibrated with 
PIA Version 9.0 / Cal G Version 6.0. For the conversion to 
an absolute flux level, measurements of the ISOPHOT internal 
Fine Calibration Source (FCS) obtained at the beginning and 
end of each scan in each filter were used. The signals from the 
two accompanying FCS measurements were linearly interpo¬ 
lated in time to calibrate each sky position of each scan in each 
filter separately. Because the ISOPHOT calibration has been 


' The ISOPHOT data presented in this paper were reduced using 
PIA, which is a joint development hy the ESA Astrophysics Division 
and the ISOPHOT Consortium. The ISOPHOT Consortium is led hy 
the Max-Planck-Institut fiir Astronomie, Heidelberg. 
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Diffuse / Localised Cirrus 



Diffuse / Localized Cirrus 


Diffuse / Localized Cirrus, Zodical Light 




Diffuse / Localized Cirrus 




Diffuse / Localized Cirrus, Zodical Light 


Diffuse + Localized Cirrus + Zodical Light 




Diffuse / Localized CirruS, Zodical Light 



Radial Distance 


Radial Distance 


Radial Distance 


Fig. 1. One-dimensional model of a FIR scan measurement 
containing cirrus only. A diffuse (straight line) and a localized 
cirrus component (gaussian) at 120 pm (top left) and 180 pm 
(top right) have similar temperatures, hence the summed inten¬ 
sities at 120 pm (bottom left, continuous) and 180 pm (bottom 
left, dashed) have a constant ratio (bottom right). 


updated since the publication of the Coma results (Stickel et al. 
1998) it was re-analyzed to get a homogeneous set of data. 

A comparison of the various signal derivation methods 
showed that the most robust results as judged by the scatter 
in the scan signals and Ii 20 Aim/Ii 80 Aim surface brightness ra¬ 
tios were obtained with the trimmed pairwise distribution us¬ 
ing the myriad value as central estimator, while the median and 
mean values of the trimmed distribution showed a larger scat¬ 
ter. Therefore, only the calibrated data derived from the pair¬ 
wise myriad method will be considered further. However, the 
gross features in the Ii 20 Mm/Ii 80 Atm surface brightness ratios 
described below were present with decreasing confidence in 
the Ii 20 Atm/Ii 80 Atm profile of all different methods used. 


5. Data Analysis 

To search for systematic trends in the FIR color profiles across 
the clusters, the Ii 20 Mm/Ii 80 Atm ratios for each position angle 
were computed for each detector pixel separately and subse¬ 
quently averaged. If necessary, at most a single outlier at each 
sky position was rejected. To further reduce the scatter in the 
Ii 20 Atm/Ii 80 Atm ratios, peaks from point sources were removed 
by linearly interpolating between the two adjacent sky posi¬ 
tions. The overall flux ratio distribution was derived by averag¬ 
ing the individual ratios of all pixels and both PAs, applying a 
global rescaling or outlier rejection if necessary. 

Care has to be taken if in addition to the diffuse large scale 
FIR emission from the galaxy a localized cirrus structure hav¬ 
ing a similar temperature is crossed. Even the small contribu¬ 
tion from the at FIR wavelengths weak zodiacal light will break 
the degeneracy of the diffuse and localized cirrus component in 
the Ii 20 /im/Ii 80 Aim T^tio and mimic an additional extended FIR 
emitting component. 


Fig. 2. One-dimensional model of a FIR scan measurement 
with cirrus and zodiacal light. In addition to the diffuse (con¬ 
tinuous straight line) and a localized cirrus component (gaus¬ 
sian), a zodiacal light component (dashed line) with a much 
higher temperature is present at 120 fim (top left) and 180 fim 
(top right). Although the summed intensities (bottom left) look 
similar to the case without zodiacal light, their ratio (bottom 
right) does reflect the profile of the localized cirrus structure as 
a dip. 


Specifically, if the 120 pm and 180 pm surface bright¬ 
nesses of the diffuse cirrus component are related by 
Ii 20 = X Ii 80 ’ brightnesses of a localized 

cirrus structure with a similar temperature are given by 
Il° 2^0 = R"" X I‘i° 80 . the Ii 20 ;.m/Ii 80 Mm ratio of galactic 
cirrus alone is given by 


jdif 

'120 


jloc 

-^120 


jdif I Jloc 

'ISO "T '180 




( 1 ) 


For similar temperatures of both components, this ratio is 
nearly constant, and almost independent of the brightness pro¬ 
files of both cirrus components. For an additional zodiacal light 
component the brightnesses at both wavelengths are related by 
If°di = X Ififj', where ^ j^cirr because of the 

different temperatures of galactic cirrus and zodiacal light. 

The total observed Ii 20 Atm/Ii 80 Aim ratio is then given by 

jdif 


jaij I jloc 
^120 -^120 


Tzodi 
■ -^120 


jdif I jloc 
-^180 -^1 


^180 


Tzodi 
■ -^180 


= i?" 


^j^zodt _ j^cirr^ 


Tzodi 

-^180 


jdif 

^180 


jloc 

-^180 


Tzodi 
■ -^180 


( 2 ) 


Although the brightnesses of the diffuse cirrus component 
and the zodiacal light vary only slightly across the cluster, the 
overall ratio is not constant due to the varying brightness of 
the localized cirrus structure I'l'^OAtm. ^rid directly reflects the 
brightness profile of this component. 

This behavior can be exemplified by a one-dimensional 
model, where the diffuse cirrus component I*^*^ might be ap¬ 
proximated by a slanted line with a constant ratio 1^*20 / IigO’ 
while the localized cirrus structure might be modeled for sim- 
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plicity by a gaussian with a constant I*° 2 o ! I'lso ratio equal to 
that of the diffuse component. The zodiacal light component 
can also be described by a slanted line to account for a gradient 
in the zodiacal light distribution along the scan path. 

While the model profile of (Eq. is constant (Fig. 
1^, the zodiacal light component in makes 

the presence of a localized cirrus component obvious, leading 
to a characteristic change in the Ii 20 /im/Ii 80 /im ratio along the 
model scan, which actually reflects the profile of the localized 
cirrus component (Fig. 

On the other hand, if no localized cirrus struc¬ 
ture but instead an ICD component described by 
Ii 2 o)tm ~ X is present with a temperature 

different from the galactic cirrus, i.e. ^ then the 

ratio 



[R^^^ - X 


tICD 

-^180 

jdif , jiCD 
-^180 -^180 


(3) 


is already varying across the cluster due to the non-constant 
ICD component, reflecting the ICD brightness profile. Adding 
the zodiacal light only changes the level of variation in the ob¬ 
served overall Ii 20 Aim/Ii 80 /im ratio, but not the presence of a 
dip or bump originating in the ICD component. 

This behavior provides the possibility to distinguish in the 
zodiacal light subtracted Ii 20 Atm/Ii 80 Atm ratio an additional 
ICD component with a surface brightness varying on scales 
of a clusters diameter from a localized cirrus component. If 
a dip or bump is still present in the Ii 20 /im/Ii 80 /im ratio after 
the subtraction of the zodiacal light component, it can be at¬ 
tributed to an additional FIR component outside the Galaxy 
and most likely identified with the thermal FIR emission of 
ICD. On the other hand, a dip or bump which disappears in the 
Ii 20 Aim/Ii 80 Aim ratio with the subtraction of the zodiacal light 
component is most likely due to a localized cirrus structure. 
Only if the ICD properties (temperature, grain composition, 
etc.) closely resemble that of galactic cirrus, a characteristic 
change of the Ii 20 Atm/Ii 80 Atm ratio from ICD will also disap¬ 
pear with the subtraction of the zodiacal light component, and 
there is then no possibility to separate the two components with 
the measurements at two FIR wavelengths. 


6. Results 


6.1. Abell 262 

Abell 262 is a spiral rich cluster with an E/D type central galaxy 
(NGC 708) showing dust lanes and diffuse gas with line emis¬ 
sion. The cluster X-ray temperature is rather low, and there is 
independent evidence for a central cooling flow (Tab. |]). Gio- 
vanelli & Haynes (1985) found that a significant fraction of 
the galaxies inside one Abell radius are deficient in HI, while 


Bravo-Alfaro et al. (1997) have shown that a number of these 


galaxies have an asymmetric HI distribution. Both results are 
most likely the results of the interaction of the cluster galaxies 
with the ICM. 



Distance [arcmin] 



Distance [arcmin] 

Fig. 3. The observed brightness distributions of the four C200 
detector pixels at 120 ^m for Abell 262 along PA 45° (top) 
and PA 135° (bottom). South lies to the left at negative dis¬ 
tances, north to the right. The brightness level is correct only 
for the lowest data stream. Eor clarity, the other three pixel data 
streams are offset arbitrarily. The 180/rm brightness distribu¬ 
tions (not shown) are quite similar. 


The 100 ^m IRAS HIRES image shows an elongated 
cloud-like patch of EIR emission, roughly oriented at a PA 
sa 50°, overlaid by filamentary structures. One ISOPHOT scan 
was executed parallel to this general elongation, while the other 
one was perpendicular. The extended 100 pm emission region 
lies at the edge of a much larger cirrus structure, and David 
et al. (1996) argued that the excess X-ray absorber required 
by the ROSAT PSPC observations is actually galactic cirrus 
rather than an absorber physically associated with the cluster. 
The inferred column density of galactic HI is about twice the 
value usually used for the analysis of X-ray data (Stark et al. 
1992), indicating significant fluctuations in the galactic HI col¬ 
umn density on angular scales < 2°. 

However, from emission line ratios ^u (1992 ) found evi¬ 
dence for an excess reddening of Eb-v ~ 0. 20 mag, indicating 
a cluster internal absorber containing dust. Wise et al. (1993 ) 
attempted to subtract the galactic foreground cirrus from IRAS 
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Distance [arcmin] 


Fig. 4. The raw surface brightness ratios Ii 20 /im/Ii 80 /im, aver¬ 
aged over all four detector pixels, along PA 45° (asterisks) and 
PA 135° (squares) as a function of distance from the cluster 
center of Abell 262. South lies to the left at negative distances, 
north to the right. Both PAs show a marked south-north asym¬ 
metry and a broad depression closely aligned with the cluster 
center. 



Fig. 5. The surface brightness ratios Ii 20 Atm/Ii 80 /im after sub¬ 
traction of the Zodiacal light, averaged over all four detector 
pixels, along PA 45° (asterisks) and PA 135° (squares) as a 
function of distance from the cluster center of Abell 262. South 
lies to the left at negative distances, north to the right. 


ISSA plates and found evidence for residual diffuse extended 
emission at 60 p,m and 100 pm, which was attributed to intra¬ 
cluster dust, too. 

Both ISOPHOT scans show for both position angles the 
broad bump of extended FIR emission, and confined to one 
sky position, the presence of an unresolved central point source 
(Fig. It can not be decided whether this FIR point source is 
due to the central cD g alaxy, because a spiral galaxy is pro¬ 
jected onto its ou tskirts (j^anti et al. 1982 ). This confusion was 
already noted by Bregman et al. (1990|). Several other compact 


Fig. 6. The overall zodiacal-light subtracted surface brightness 
ratio Ii 20 /im/Ii 80 /im for Abell 262, averaged over both position 
angles and all detector pixels. 


sources appear only in the data streams of individual pixels, 
indicating non-central crossings. 

The raw Ii 20 /xm/Ii 80 Aim surface brightness ratios including 
zodiacal light for both position angles, separately averaged over 
the four detector pixels, agree remarkably well (Fig. Q), even in 
the marked north - south asymmetry, and show a very broad 
minimum somewhat off-center from the nominal cluster centre 
position. Remarkably, the broad depression completely disap¬ 
pears in the zodiacal-light subtracted surface brightness ratios 
along both PAs (Fig.^), leaving only an almost linearly increas¬ 
ing overall Ii 20 A.m/Ii 80 Mm ratio (Fig. H). 


6.2. Abell 400 


A dynamical analysis of Abell 400 was carried out by Beers et 
al. (1992), who concluded that this spiral rich cluster actually 
consists of two bound sub-clusters currently undergoing merg¬ 
ing. The irregular X-ray morphology ( Beers et al. 1992| ; Buote 
& Tsai 1996) and the absence of a cooling flow (White et al. 


1997; Allen & Fabian 1997) also indicates a young dynami¬ 
cal cluster age. The overall X-ray temperature of w 2.5 keV is 
rather low (David et al. 1993; [White et al. 1997| ). 

Although Abell 400 lies at a much higher galactic latitude 
than e.g. Abell 262, both ISOPHOT scans (Fig. ^ show nev¬ 
ertheless a much higher absolute level of the 120/rm surface 
brightness. The 100 pm IRAS HIRES image shows quite a 
structured region with large patches of cirrus foreground to¬ 
wards Abell 400. Two very bright point sources were detected 
with one detector pixel only along PA 40° and PA 130 °, which 
have been identified with PGC 11141 (IRAS F02542H-0600) 
and UGC 2444 (IRAS 02558H-0606), respectively. Both scans 
also indicate in the data streams of several pixels the presence 
of a much weaker compact source confined to the sky position 
centered on the cluster. 

The Ii 20 Atm/Ii 80 A»m raw surface brightness ratios (Fig. |) 
along the two PAs are quite structured on angular scales of 
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Distance [arcmin] 


Fig. 7. The observed brightness distributions of the four C200 
detector pixels at 120/rm for Abell 400 along PA 40° (top) 
and PA 130° (bottom). South lies to the left at negative dis¬ 
tances, north to the right. The brightness level is correct only 
for the lowest data stream. For clarity, the other three pixel data 
streams are offset arbitrarily. The 180 pm brightness distribu¬ 
tions (not shown) are quite similar. 


10' - 20' without a good overall agreement of the two pro¬ 
files. However, the ratio at the scan crossing is in almost per¬ 
fect agreement, thereby providing an independent check of the 
flux calibration. After subtraction of the zodiacal light, the 
Ii 20 /xm/Ii 80 Aim surface brightness ratios (Fig. ||) of the two PAs 
are in much better agreement, showing a flatter part towards 
south and a gradual steepening towards north. This overall be¬ 
havior is quite easily seen in the overall Ii 20 Aim/Ii 80 Aim surface 
brightness ratio (Fig. p^. 


6.3. Abell 496 


The relaxed single symmetric X-ray morphology (Buote & 
Tsai 1996), the significant cooling flow (Tab. |^, and the 
double- nucleus centr al cD galaxy with very small velocity dis¬ 
persion ( Tonry 1985 ) indicates a dynamically old, relaxed clus¬ 
ter. There is evidence for an X-ray absorption column density 



Fig. 8. The raw surface brightness ratios Ii 20 /xm/Ii 80 Mm aver¬ 
aged over all four detector pixels, along PA 40° (asterisks) and 
PA 130° (squares) as a function of distance from the center of 
Abell 400. The profiles along the two PAs are rather dissimilar 
with structures on angular scales of 10' - 20'. South lies to the 
left at negative distances, north to the right. 



Distance [arcmin] 

Fig. 9. The surface brightness ratios Ii 20 /im/Ii 80 Aim after sub¬ 
traction of the Zodiacal light, averaged over all four detector 
pixels, along PA 40° (asterisks) and PA 130° (squares) as a 
function of distance from the center of Abell 400. South lies 
to the left at negative distances, north to the right. The overall 
shapes of both PAs show a flatter region towards south, and a 
steepening towards north. 


above the galactic value across the cluster (White et al. 1994; 
|VlacKenzie et al. 1996 ; Allen et al. 2001), consistent with the 
excess reddening of Eb-v ~ 0.20 mag derived from emission 
line ratios (Hu 1992). 

The ISOPHOT scans along both PAs (Fig. 0 ) are devoid 
of strong point sources. In the data streams of two pixels along 
PA 180°, there is an indication for a weak compact source con¬ 
fined to a single sky position. A broad hump offset to the east 
from the cluster center is visible along PA 90°, and a similar 
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Fig. 10. The overall zodiacal-light subtracted surface bright¬ 
ness ratio Ii 20 Aim/Ii 80 /im fof Abell 400, averaged over both 
position angles and all detector pixels. 




Distance [arcmin] 


Fig. 11. The observed brightness distributions of the four C200 
detector pixels at 120/rm for Abell 496 along PA 90° (top) 
and PA 180° (bottom). The brightness level is correct only for 
the lowest data stream. For clarity, the other three pixel data 
streams are offset arbitrarily. The 180/rm brightness distribu¬ 
tions (not shown) are quite similar. 



Distance [arcmin] 


Fig. 12. The raw surface brightness ratios Ii 20 /im/Ii 80 Aim, av¬ 
eraged over all four detector pixels, along PA 90° (asterisks, 
from east towards west) and PA 180° (squares, from south to¬ 
wards north) as a function of distance from the center of Abell 
496. The prohle is nearly constant along PA 90°, while PA 180° 
shows an almost linear increase. 



Distance [arcmin] 

Fig. 13. The surface brightness ratios Ii 20 Atm/Ii 80 Aim after sub¬ 
traction of the Zodiacal light, averaged over all four detector 
pixels, along PA 90° (asterisks, from east towards west) and 
PA 180° (squares, from south towards north) as a function of 
distance from the center of Abell 496. 


but weaker bump appears to be superposed on the generally 
decreasing surface brightness profiles along PA 180°. 

The two raw Ii 20 Aim/Ii 80 Aim surface brightness ratios in¬ 
cluding the zodiacal light (Fig. along PA 180° and PA 90° 
are strikingly different. The former is almost monotonically in¬ 
creasing from south towards north across the cluster, while the 
latter varies much less across the cluster. The color prohles af¬ 
ter subtraction of the zodiacal light (Fig. do not change 
very much, still showing a steady increase from south to north 
and a rather flat but wiggling behavior along PA 90°. Neverthe¬ 
less, both prohles were combined to check for the presence of 
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Distance [arcmin] 

Fig. 14. The overall zodiacal-light subtracted surface bright¬ 
ness ratio Ii 20 /im/Ii 80 /im for Abell 496, averaged over all four 
detector pixels. 


structure on scales of 10' - 20'. The resulting profile (Fig. 
however, shows only a gentle rise. 


6.4. Abell 1656 ( Coma) 

There is now overwhelming evidence for recent or ongoing 
merging processes in the Coma cluster from optical and X-ray 
data ( Colless & Dunn 1996|; |Vikhlinin et al. 1997|; B urns et al. 
1994b; [Neumann et al. 2001|). The early evid ence for intraclus¬ 
ter light ( Welch & Sastry 1971 ; Mattila 1977 ) has recently been 
confirmed with the detection of additional low surface bright¬ 
ness struct ures interpreted as tidal debris from recent galaxy in - 
teractions ( [Gregg & West 19% ; Trentham & Mobasher 1998 ). 
With an X-ray temperature of w 8 keV, Abell 1656 is the hottest 
of the observed clusters, and in fact lies near the upper end of 
the X-ray luminosity - temperature correlation (Tucker et al. 
1998). 

Except for a slight decrease in surface brightness level due 
to the updated ISOPHOT calibration, the scans along PA 82° 
and PA 36° (Fi g, jf^ appear almos t unchanged compared to 
those shown in Stickel et al. (1998). However, due to the im¬ 


proved data processing, there is now a much closer agreement 
between the Ii:^m/Ii 80 Atm surface brightness ratios along the 
two PAs (Fig. |I^, showing significant deviations only at the 
outer ends of the scans. An increased Ii 20 Aim/Ii 80 /im ratio 
within « 10' from the cluster center is apparent, which essen¬ 
tially remains unchanged even after the zodiacal light has been 
subtracted (Fig. 0. The overall Ii 20 /im/Ii 80 um ratio averaged 
over both PAs and all detector pixels (Fig. shows the bump 
quite clearly, which can be quite well represented by a gaussian 
with a peak 5% above a somewhat slanted background. 


6.5. Abell 2670 


The X- ray morphology of Abell 2670 has a double peak in the 
center (Buote & Tsai 1996|), where the secondary peak is due to 



Distance [arcmin] 



Distance [arcmin] 


Fig. 15. The observed brightness distributions of the four C200 
detector pixels at 120 p,m for Abell 1656 along PA 82° (top) 
and PA 36° (bottom). The brightness level is correct only for 
the lowest data stream. For clarity, the other three pixel data 
streams are offset arbitrarily. The 180/im brightness distribu¬ 
tions (not shown) are quite similar. 


several point sources, and a pronounced asymmetry at low flux 
levels almost along the north-south direction (Hobbs & Will- 
more 1997). The ISOPHOT scan position angles were chosen 
to cover both features. 

There are conflicting results regarding the dynamical stage 
of Abell 2670. [Bird (1994 , ^scalera etal. (1994 , and Hobbs & 
Willmore (1997) argue in favour of significant substructure as 
a result of merging, while Sharpless et al. (1988) and Girardi et 
al. (1995) pointed out the contamination by a foreground clus¬ 
ter. This controversial picture is reinforced by the analysis of 
X-ray data, where conflicting results about a possible cooling 
flow have been found (see Tab. ||). 

Extended FIR emission was described by Wise et al. 
(1993), but judged as uncertain due to the relatively large clus¬ 
ter redshift and small angular extent, and also by Cox et al. 
(1995). This diffuse emission is easily seen on the 100 p,m 
IRAS HiRes image. However, rather than being centered on 
the cluster, it covers only the western part of the cluster with 
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Fig. 16. The raw surface brightness ratios Ii 20 Atm/Ii 80 Atmc av¬ 
eraged over all four detector pixels, along PA 82° (asterisks) 
and PA 36° (squares) as a function of distance from the center 
of Abell 1656. Both distributions show a broad bump closely 
aligned with the cluster center, indicating a cold component 
with Il20Aim > Il80/im- 



Distance [arcmin] 


Fig. 17. The surface brightness ratios Ii 20 /im/Ii 80 /im after sub¬ 
traction of the Zodiacal light, averaged over all four detector 
pixels, along PA 82° (asterisks) and PA 36° (squares) as a func¬ 
tion of distance from the center of Abell 1656. The profiles are 
almost unchanged compared to Fig. |^. 


a north-south elongated patch of low level flux, overlaid with 
some slightly resolved brighter structures. The central compo¬ 
nent of these structures is seen in the middle of the overall 
rather flat surface brightness distributions along both PA 84° 
and PA 19° (Fig. |i|). 

The Ii 20 /im/Ii 80 /im surface brightness ratios (Fig. from 
both PAs show a general similarity, a narrow central dip closely 
aligned with the cluster centre, increasing to approximately 
the same level at the outer edges of the scans. After subtrac¬ 
tion of the zodiacal light, the central dip vanishes in both in¬ 
dividual scans (Fig. ^l|), resulting in an almost flat overall 


Fig. 18. The overall zodiacal-light subtracted surface bright¬ 
ness ratio Ii 20 /im/Ii 80 /jm for Abell 1656 averaged over both 
position angles and all detector pixels. 



Distance [arcmin] 



Distance [arcmin] 


Fig. 19. The observed brightness distributions of the four C200 
detector pixels at 120p,m for Abell 2670 along PA 84° (top) 
and PA 19° (bottom). The brightness level is correct only for 
the lowest data stream. For clarity, the other three pixel data 
streams are offset arbitrarily. The 180/tm brightness distribu¬ 
tions (not shown) are quite similar. 
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Fig. 20. The raw surface brightness ratios Ii 20 Atm/Ii 80 Atm, av¬ 
eraged over all four detector pixels, along PA 84° (asterisks) 
and PA 19° (squares) as a function of distance from the center 
of Abell 2670. 



Distance [arcmin] 

Fig. 22. The overall zodiacal-light subtracted surface bright¬ 
ness ratio Ii 20 /im/Ii 80 /jm for Abell 2670 averaged over both 
position angles and all detector pixels. 



Distance [arcmin] 

Fig. 21. The surface brightness ratios Ii 20 /im/Ii 80 /im after sub¬ 
traction of the Zodiacal light, averaged over all four detector 
pixels, along PA 84° (asterisks) and PA 19° (squares) as a func¬ 
tion of distance from the center of Abell 2670. 


Ii 20 Atm/Ii 80 Atm ratio across the cluster with only small scale 
low-level variations (Fig. 


6.6. Abell 4038 (Klemola 44) 


The X-ray morphology of Abell 4038 is similar to that of 
the Coma cluster, being essentially elliptical with a faint ex¬ 
tended source representing a group of galaxies within the clus¬ 
ter (Buote & Tsai 1996). The X-ray data indicate the presence 
of a cooling flow (see Tab. The overall galaxy distribution 
is elongated along PA 120°, as is the major axis of the central 
cD galaxy. (Green et al. 199C). One of the ISOPHOT scans was 
aligned with this direction while the other one was perpendic¬ 
ular. 


The overall surface brightness distributions along both PA 
120° and PA 30° (Fig. ^ are rather flat. Only along PA 120° 
compact sources show up in the data streams of single detector 
pixels. The brightest source detected along PA 120° is actually 
a background galaxy, identified with IRAS 23440-2816 (z = 
0.0495). 

The Ii20/xm/Ii80/xm surface brightness ratios (Fig. 
along the two PAs agree well, being generally flat with 
only small-scale structures. The same behavior is seen after 
zodiacal-light subtraction (Fig. and eventually in the over¬ 
all Ii 20 Aim/Ii 80 Mm profile (Fig. 1^. 


7. Discussion 


The striking result of the analysis of the Ii 20 Aim/Ii 80 /xm 
zodiacal-light subtracted surface brightness ratios is the pres¬ 
ence of an easily noticeable bump in only one cluster, namely 
Abell 1656 (Coma), the cluster with the highest X-ray temper¬ 
ature and no indication of a cooling flow. 

Although the raw Ii 20 /im/Ii 80 Aim profiles showed struc¬ 
tures on scales of 10' - 20' also for Abell 262 and Abell 2670, 
these structures disappeared with the subtraction of the zodia¬ 
cal light component. This indicates that an additional FIR emit¬ 
ting component outside the Galaxy is present in the direction of 
Abell 1656 while galactic cirrus foreground emission close to 
the line-of-sight is responsible for the structures in the color 
profiles of Abell 262 and Abell 2670. This is in clear con¬ 


trast to the findings of Wise et al. (1993), who ascribed resid¬ 


ual FIR emission towards Abell 262 and Abell 2670 found on 
IRAS ISSA plates to ICD. However, in these two as well as the 
other three observed clusters, the presence of ICD cannot be 
excluded completely, because it might be present at very low 
levels or have properties similar to galactic cirrus. In the for¬ 
mer case, the signature is easily lost in the noise, while in the 
latter case there is no way of discriminating it against the fore- 
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Distance [arcmin] 



Distance [arcmin] 

Fig. 23. The observed brightness distributions of the four C200 
detector pixels at 120 /im for Abell 4038 along PA 120° (top) 
and PA 30° (bottom). The brightness level is correct only for 
the lowest data stream. For clarity, the other three pixel data 
streams are offset arbitrarily. The 180/im brightness distribu¬ 
tions (not shown) are quite similar. 

In the scan along PA 30° (bottom), one pixel data stream shows 
a marked jump in surface brightness near the cluster center, not 
seen in any other pixel or wavelength. This is most likely an 
abrupt change in pixel sensitivity not accounted for by the cal¬ 
ibration. This pixel data stream has therefore been ignored in 
the subsequent analysis. 



Fig. 24. The raw surface brightness ratios Ii 20 /im/Ii 80 Aim, av¬ 
eraged over all four detector pixels, along PA 120° (asterisks, 
average of four pixels) and PA 30° (squares, average of three 
pixels) as a function of distance from the center of Abell 4038. 



Fig. 25. The surface brightness ratios Ii 20 Atm/Ii 80 Atm after sub¬ 
traction of the Zodiacal light, averaged over all four detector 
pixels, along PA 120° (asterisks, average of four pixels) and 
PA 30° (squares, average of three pixels) as a function of dis¬ 
tance from the center of Abell 4038. 


ground galactic dust even with photometric measurements at 
additional FIR wavelengths. 

The bump in the Ii 20 /im/Ii 80 /im ratio in Abell 1656 (Coma) 
might be caused by the integrated FIR emission of the un¬ 
resolved cluster galaxies, if their ISM had different tempera¬ 
tures and/or dust properties than the foreground galactic cir¬ 


1656 were rather special, the other clusters should also show a 
characteristic signature in their Ii 20 /im/Ii 80 Aim profiles, which, 
however, is not observed. Moreover, the FIR color temperature 
determination of a large sample of norm al galaxies selected 
from the ISOPHOT Serendipity Survey (Stickel et al. 200C) 


showed the ubiquitous presence of dust in the ISM with tem¬ 
peratures of T 17 K similar to the Galaxy. 

In principle, the Sunyaev-Zeldovich effect might be another 
cluster intrinsic cause for extended excess emission at FIR and 
sub-mm wavelengths. However, even for very hot clusters, the 
effect at the peak wavelength of w 800/im is small (e.g. Ko¬ 
matsu et al. 1999), and vanishes shortward of « 300/im (Chase 


rus (Quillen et al. 1999). Unless the galaxy content of Abell 


et al. 1987; Rephaeli 1995), making it very unlikely as a cause 
for the observed systematic changes in the Ii 20 Atm/Ii 80 Aim pro¬ 
files. 

The most likely interpretation of the observed spatially ex¬ 
tended bump in the Ii 20 Aim/Ii 80 Aim ratio of Abell 1656 is there¬ 
fore the presence of intracluster dust mixed with the hot X- 
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Distance [arcmin] 


Fig. 26. The overall zodiacal-light subtracted surface bright¬ 
ness ratio Ii 20 /im/Ii 80 Aim for Abell 4038 averaged over both 
position angles and all detector pixels. 


ray emitting electron gas. The excess surface brightness in¬ 
dicates that the intracluster dust has properties different from 
the foreground galactic cirrus such as temperature, grain sizes, 
grain composition, or emissivity, otherwise no signature in the 
Ii 20 Aim/Ii 80 Aim ratios would be discernible. 

The observed overall zodiacal light subtracted 
Ii 20 Mm/Ii 80 Mm r^tio Can be written as 

Tcirr _i_ tICD 
yobs ^120 -^120 


j^OOS _ 


(4) 


Tcirr I jICD ’ 

-^180 “'"'^180 

where Ij 2 q and are the surface brightnesses of the fore¬ 
ground galactic cirrus and the excess emission of the cluster 
at 120 ^m, respectively, and IlfgQ and the corresponding 
brightnesses at 180/rm. To quantify the bump in the overall 
zodiacal-light subtracted Ii 20 /im/Ii 80 Aim surface brightness ra¬ 
tio of Abell 1656, it was modelled as the sum of a gaussian, 
presumed to result from the cluster dust emission and a 
first order polynomial from the large scale cirrus foreground 
fcirr rpjjg cin-us foreground can then be interpolated at the po¬ 
sition of the cluster centre, giving the central ISr2o/I585 ratio 


Tcirr 

ycirr -^120 


Tcirr 

■^180 


j^cirr _ 


(5) 


of the cirrus alone. 

For the case of > R""", 
can be re-arranged to give 


the two relations (Eqs. 


tICD 

'120 


= r 


180 


(^ob. _ 


tICD 

r\ _|_ jyobs -^180 
/ ' Tcirr 

-^180 


(6) 


Under the assumption that the cluster dust emission is weak 
compared to the cirrus foreground, i.e. I^gJ* <C Iig'(5, it can fur¬ 
ther be rearranged to 

(7) 


(TCirr I tICD\ ( ryobs rycirr\ 
- UlSO + ^180 ) ^ ^ )■ 


tICD 
^120 

This gives the cluster excess surface brightness 1^2? in terms of 
the observed overall Ii 20 /im/Ii 80 /im ratio R°'^®, the interpolated 
Ii 20 Aim/Ii 80 Aim ratio R'^*’''' of the cirrus, and the observed total 
surface brightness at 180 p,m. 


For Abell 1656, R°'"" - R"" = 0.07, with a characteristic 
e-folding angular scale of cticd ~ 5', and I^g*^ + l{gj^ = 3.1 
MJy/sr, resulting in a cluster excess peak surface brightness 
Ii 2 ? « 0.2 MJy/sr. Assuming that the fitted one-dimensional 
gaussian is representative of an underlying circular symmet¬ 
ric two-dimensional excess flux distribution, the integrated (to¬ 
tal) excess fluxes computed from Ficd = 2 ttctjcd the 
volume of the two-dimensional gaussian with a peak height 
equal to the derived excess surface brightness and a width 
equal to the fitted one-dimensional gaussian, is Fjcd ~ 
2.8 Jy. This corresponds to a monochromatic FIR luminosity 
Licd = 47rD^:/FicD ofLicD ~ 8.6 x lO'^^erg/s, whereD is 


the cluster luminosity distance. 

The absence of a clear sign of cluster excess emission in 
the individual scans of Abell 1656 at both wavelengths gives 
strict upper limits for the ICD excess surface brightness of 
1^20* < 0.4MJy/sr. Under the assumption that the profile of 
the ICD excess emission at both wavelengths can be approxi¬ 
mated by a gaussian and the cirrus foreground surface bright¬ 
ness varies linearly across the cluster, this absolute upper limit 
can be transformed to a lower limit of for the peak surface 
brightness ratio ij^f /Iis? of foo excess emission. This is done 
by fitting the ratio of the presumed surface brightness profiles 
(for each wavelength a gaussian added to a slanted background) 
to the overall Ii 20 Atm/Ii 80 Aim ratio, subject to the constraint not 
to exceed the upper limit, yielding 1^2?/Iis? ~ 2. This lower 
limit in the lli^^/ll^? ratio transforms into a lower limit in 
the dust color temperature of T ^ 30 K, assuming a blackbody 
(Planck) function modified by a frequency dependent emissiv¬ 
ity Fj^ = z/^B^(Td) with exponent (3 = 2. 

An estimation of the intracluster dust mass can be derived 
by using 

Md=D‘^F, (8) 

( Hildebrand 1983 ; Draine 1990 ), where F;^ is the total excess 
flux density at 120/im, To the dust color temperature, and kx 
the dust opacity. With a dust opacity of 3 m^/kg, which lies in 
the middle of the range of currently considered values (Draine 
1990), a dust temperature of Td w 30 K and the excess flux of 
Ficd = 2.8 Jy a dust mass of Md ~ 10^ Mq is derived. This 
dust mass is only exemplary, since it depends strongly on the 
assumed dust color temperature and the dust opacity. 

It should be noted that a dust opacity, which already in¬ 
cluded the conversion to gas mass with a constant galactic gas- 
to-dust ratio, had been used in the first publication of the Coma 
observations (Stickel et al. 1998). The newly derived dust mass 
for Abell 1656 is therefore considerably lower than that given 
before. Consequently, the gas-to-dust ratio using a gas mass of 
Mgas ~ 10^^ Mq is extremely high (« 10®). Furthermore, the 
low inferred dust ma ss leads under standard assumptions about 
the dust properties ( Whittek 1992 ) to a negligible amount of 
visual extinction (Ay “C O.lm ag), in contras t to the reported 
optical extinction of « 0.3 mag ( Zwicky 1962 ; Karachentsev & 
Lipovet skii 1969), and in a ccord with the analysis of the IRAS 
data by Dwek et al. (1990). 
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8. Conclusions 


Spatially extended FIR excess emission from dust distributed in 
the intracluster medium has been searched for in a small sample 
of six Abell clusters covering a wide variety in cluster proper¬ 
ties, such as X-ray morphology and temperature, cluster age, 
richness, and distance, by using ISOPHOT scan measurements 
at 120/im and 180/rm wavelength. After subtraction of the zo¬ 
diacal light, the characteristic signature of a dip or bump on 
angular scales of 10' - 20' is only seen in Abell 1656 (Coma), 
a hot, dynamically young cluster with optical and X-ray signs 
of interactions with smaller galaxy groups. 

The evidence from IRAS ISSA plates for FIR emission 
from intracluster dust in Abell 262 and Abell 2670 (Wise et 
al. 1993) could not be confirmed. Although the raw (including 
zodiacal light) Ii 20 /im/Ii 80 /im profiles did show a dip, it dis¬ 
appeared after the zodiacal light had been subtracted. This be¬ 
havior is indicative of a localized patch of galactic foreground 
cirrus, which happens to lie along the line of sight to these clus¬ 
ters. This strikingly underlines the necessity of interpreting FIR 
measurements, even flux ratios at several wavelengths, care¬ 
fully by taking into consideration the ever-present cirrus struc¬ 
tures in conjunction with the zodiacal light, which can produce 
a quite counterintuitive behavior. 

This result does seem to support the simple picture put for¬ 
ward in IStickel et al. (1998 ), that only young clusters currently 
undergoing merging show the signature of intracluster dust, 
which might indicate that intracluster dust is primarily brought 
into the cluster from the outside, either by merging with exter¬ 
nal galax y groups or possibl y by steady infall of intercluster 
material ( j^opescu et al. 20()C ). 

The FIR observations of the six clusters, particularly those 
with a strong cooling flow inferred from X-ray observations, 
do not seem to support the picture of cold clouds dropping out 
of the cooling flow and shielding intracluster dust from sput¬ 
tering, or even dust originating in the dense centers of cooling 
flow clusters. The rather low inferred dust mass in Abell 1656 
and the non-detection of the characteristic signature of ICD in 
the other five clusters can be taken as evidence that the X-ray 
excess absorption seen in cooling flow clusters is likely not due 
to the presence of intracluster dust. Particularly, no excess X- 
ray absorption is required for Abell 1656 (Allen et al. 2001), 
the only cluster where evidence for extended FIR emission has 
been found in this study. The absence of a signature for intra¬ 
cluster dust in the clusters with strong indication for cooling 
flows is in agreement with several other unsuccessful searches 
for the final stages of the cooled material, and adds additional 
doubt on the correctness of material condensations associated 
with the cooling flow model. 

Finally, intracluster dust appears not be responsible for the 
the observed dimming of the high - redshift supernovae, and the 
contribution to the FIR extragalactic background light is likely 
also vanishingly small. 
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